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ABSTRACT 

This  is  a  continuation  of  MRC  Technical  Summary  Report 
#997.  The  connection  between  monosplines  and  quadrature 
formulae  is  used  to  determine  the  best  quadrature  formula  for 
J  f(z)dz  integrated  along  the  unit  circle.  The  quadrature 
formula  exact  for  all  spline  functions  is  also  derived. 


ON  POLYNOMIAL  SPLINE  FUNCTIONS  ON  THE  CIRCLE  II. 
MONOSPLINES  AND  QUADRATURE  FORMULAE 

I.  J.  Schoenberg 

Introduction 

In  [4]  we  have  studied  the  interpolator/  properties  of  the  class  3^  ^ 

of  polynomial  spline  functions  on  the  circle  that  were  first  discussed  by 

Ahlberg,  Nilson  and  Walsh  (see  [1],  also  for  further  references) .  In  the 

present  second  paper  we  discuss  quadrature  formulae  of  the  form 

k-1 

(1)  f  f(z)dz  =  £c  f(J)  +  Rf  , 

U  0  } 

where  U  is  the  unit  circle  |z|  =1  described  counter-clockwise  and 
w  -  exp(  2Tri/k)  .  This  problem  requires  a  discussion  of  monosplines  and 
will  suggest  several  classes  of  spline  functions  and  monosplines  having  note¬ 
worthy  properties.  Free  use  is  made  of  the  results  of  [4],  in  particular  of 
B-splines  and  their  Fourier  series  expansions. 

1.  Mono^plines  and  quadrature  formulae.  Let  k  and  m  be  positive 
Integers.  Our  objective  is  to  derive  quadrature  formulae  of  the  form  (1)  which 
are  to  be  exact  for  polynomials  of  degree  not  exceeding  m-1 ,  i.e, 

(1.1)  Rf  =  0  if  U  tt  ,  . 

m-i 

This  property  is  equivalent  to  the  relations 

k-1 

(1.2)  fzvdz  =  0  =  Yj  c  Jv  ,  (v  =  0,1,...,  m-1)  , 

U  J=0  ] 
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and  these  show  that  the  quadrature  formula  (1)  has  the  property  (1.1)  if  and 
only  if 

k-1  .  m_1 

(1.3)  the  polynomial  2j  C.  xJ  vanishes  if  x  =  l,u,  ...,u> 

0  J 

We  conclude  the  following:  1.  If  m|k  then  the  polynomial  (1.3) 
vanishes  identically,  hence  =  0  for  all  j  .  2.  If  mck  then  the  poly¬ 

nomial  (1.3)  still  depends  on  k-m  arbitrary  linear  parameters.  Avoiding 
trivial  cases  we  shall  assume  that 


(1.4) 


1  £  m  <  k 


Expanding  the  binomials  appearing  in  (1 . 5)  we  see  that  the  condition 
(1 . 3)  is  equivalent  to  the  validity  of  the  identity  in  z 

(1.5)  £  C  (z-J)1""1  =  0  . 

0  ] 

At  this  point  we  reintroduce  the  class  S  .  of  polynomial  spline  functions  on 

m,  K 

U  of  degree  m-1  having  the  <J  as  knots,  and.  for  a  good  reason.  In  fact  the 
identity  (i .  5)  is  the  precise  condition  for  the  existence  of  a  function  S(z)c8^ 
such  that 

(1.6)  Cj  =  the  jump  of  S^m  ^(z)  at  z  =  (J  ( j  =  0,  .  .  . ,  k— 1 ) 

k-1 

Indeed,  if  S(z)  =  Q(z)«  it  ^  on  the  arc  (w  ,1)  then  the  function  on  U 

k-1  k-1 

defined  on  the  arcs  (u>  ,  1 ),  ( 1,  w) ,  . . . ,  ( w  ,1)  by  the  successive  partial 


sums  of  the  expression 


Q(  z)  +  I  C 


k-1  j.m-1 

V  (z-w') 


(m-1)! 


is,  in  view  of  (1 .  5) ,  a  spline  function  S(z)  with  the  property  (1.6) . 

The  connectior  of  the  spline  function  S(z)  with  the  quadrature  formula 
(1)  is  close  as  shown  by 

Lemma  1 .  If  we  write 
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I 


(1.7) 


K(  z)  =^T2m- S(  z) 
m  i 

and  If 

(1.8)  f(  z)  «  Cm  ( U) 
then  we  have  the  relation 

k-1  . 

(1.9)  f  f(  z)  dz  =  V  C  f(  J)  +  (-l)m  /*  K(  z)  f  m  ( z)  dz 

U  0  J  u 

Proof.  By  (1.6)  and  (1.7)  we  find  that 

(1.10)  -  Cj  =  the  Jump  of  K^m  ^(z)  at  z  -  <J 


Integrations  by  parts  show  that 


f  K(  z)  f(m)(  z)  dz  =  -/ K,(z)f(m_1)(z)dz  =  ...  =  (-l)rn'1/K(m'1)(z)f,(z)dz 
U  U  U 

a  point  beyond  which  we  can  not  continue  the  process  because  of  the  discon- 

tinuities  of  K'  '(z)  .  However,  we  may  proceed  as  follows 

J+l 

f K(z)f(m)(z)dz  =  (-l)m~l£  /  K(m'U(z)f'(z)dz 


J  — 


J+l 

K(m-1)  J+l  „(m-l).  J.,.  J.  ?  ..  .  .„(m-l).  . \ 

K  (w_  )  -  K  (w+)f (w  )  -  /  f( z)  dK  (z)j 

or  ./ 


and  using  (1. 10)  we  obtain 

f  v  /  f(^)  /  - 


/  K(z)  f  m)(t)  <.iz  =  (-1)"1'1  ^C.f(J)  +  (-  l)m/f(z)dz 

TT  >  TT 


because  d  K^m  ^(z)  =  dz  within  each  of  the  arcs. 


Using  a  term  familiar  from  the  real  axis  case  (see  e.g.  [3])  we  call  K(z) 

a  monospline  of  degree  m  and  denote  their  class  by  the  symbol  ft  .  We 

m,  k 

may  then  summarize  our  results  as  follows:  To  every  quadrature  formula  (1) 
with  the  property  (1.1)  corresponds  a  monospline  K(z)  producing  the  identity 
(1.9).  Conversely,  to  every  monospline  K(z)  corresponds  a  quadrature 
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formula  (1),  the  being  described  by  (1.6).  This  correspondence  is  one- 

to-one  up  to  an  element  of  it  .  that  we  may  add  to  K(z)  without  changing 

m-i 

our  results.  Our  discussion  of  the  quadrature  formula  (1)  will  single  out  cer¬ 
tain  spline  functions  and  monosplines  for  special  study. 

2.  The  flower-shaped  spline  functions.  Let  S(z)«  8  and  let  us 

““_J1  «  3B=a—a  ii  ii  ■  m,  k 

assume  that  it  satisfies  the  functional  equation 

(2-1)  S(zu)  =  w's(z)  ,  (|z|=l)  , 

where  r  is  an  integer,  0  ir  S  k-1  .  We  could  say  that  S(z)  is  quasiperiodic, 
or  r-quasiperiodic.  However,  in  view  of  the  rotational  symmetry  of  the  image 
of  U  by  w  =  S(z)  we  prefer  to  say  that  S(z)  is  flower- shaped  and  also  that 
S(z)  is  an  r-flower  (see  Fig.  1  in  [4]  showing  the  image  of  U  by  a  3-flower 
for  k  =  5  ,  m  =  2  )  . 

In  terms  of  the  B-spline  M  (z)  of  f  4]  we  have 

Lemma_^.  For  each  r  ,  0  ir  s  k-1  ,  there  is  up  to  a  constant  factor  a 

unique  r-flower  given  by  the  formula 

k-1 

(2.  2)  S( z)  =  Y.  JrM  (zO  . 

J=0  m 

Proof:  By  [  4  ,  Lemma  1  ]  we  may  write 

S(z)  =  £ c.  M  (zw'1)  , 

J  m 

and  (2.1)  gives  Tx.  M  (zw1  *)  =  71c.  (zu  *)  or  71C,,.M  (zw-^)  = 

LJ  J  m  LJ  j  m  LJ  Jfl  m 

Z^JmJzu,  ,  whence  =  c^o/  and  finally  =  cQ<Jr  .  Setting 
c0  =  '  we  obtain  (2.2)  .  Conversely,  it  is  easily  established  that  (2.2) 
satisfies  (2.1)  and  Lemma  2  is  established. 

The  Fourier  expansions  of  r-flowers  are  described  by 
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(2.  3) 


Lemma  3.  Jf  S(z)  is  an  r-flower  then,  up  to  a  constant  factor. 


S(z)=z  if  r  =  0, 1, . . . ,  m-1  , 


and 

(2.4) 


ks+r 


S(z)  =  2j  (ks+r)(ks+r-l).. .  (ks+r-m+1) 
s=-°° 

if  r  =  m,  m+1,  ,  k-1 

Proof:  We  recall  (4,  formula  (2.16)]  that 


oo 


(2.5) 

where 

(2.6) 


M 


Im(Z)  =  2rri  ^  bv-i  bv- 


-oo 


m+1 ' 


b 

v  v 


—  V 


if  v  #  0  ,  bQ  =  ^r 


From  (2.2)  and  (2.5) 


2.1  Z\---  b.-m*ll)=() 


k-1  .  .  . 

£  Jr  -Jv  ] 


U  U) 


or 

(2.7) 


S(  z)  =  £b  b  ,  ..  .  b 

2ni  ^  v  v/-l  v-m+1 
v 


where  the  summation  is  over  all  v  =  r  (mod  k)  ,  hence  v  =  ks  +  r  .  By  (2.6) 
we  see  that  b^  =  0  if  v  is  a  nonvanishing  multiple  of  k  .  If  0  s  r  5  m-1 
then  (2.7)  is  seen  to  reduce  to  the  single  term  corresponding  to  s  =  0  or 
v  =  r  and  this  is  indeed  proportional  to  (2.3).  Also  directly  it  is  clear  that 
and  that  z  i«  an  r-ficwer. 

If  m  %  r  §  k-1  ,  we  may  discard  on  the  right  side  of  (2.7)  the  non¬ 
vanishing  factor 

_ k.,i  ~rwi  ~r+1»  m  -r+m-1 

~  ( 1-w  )  ( 1— w  ) .  .  .  (1-w  ) 


z  €  g 


m,  k 


and  we  obtain  the  expansion  (2.  4) 
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3.  Thejnonosj^ne^Mea^t  L^-norm.  We  have  already  used  the  fact 


f’  u 

(3.  1) 


k-1 

S(z)  =  Yj  c,  M  (z  u"f 
j  m 


0 


'-eDresentf  the  most  general  element  of  8  .  .  Using  (2.  5)  we  find  its 

m,  k 


Fourier  expansion  to  be 


k-1 


S(z)  =  JL  V  b  b  ,  . . .  b  zv  Y  c 

2rri  „  v  v  — 1  v-m+1  ,un  . 


-vj 


CO 


1=0 


However 
(3  >) 


k-1 

=  Tj  ciw 

1=0  J 


-vj 


in  the  finite  Fourier  series  representation  of  an  arbitrary  periodic  sequence 
(  tj  )  of  period  k  .  This  establishes 

Lemma  4.  The  most  general  element  of  8^  ^  is  given  by 

(3.  3) 


00 


S(z)  =  —  7,  b  b  ,  . . .  b 

2rri  u  v  v-1  v- 


-00 


i  h  z 

v-m+1  v 


where  ( r\  )  is  an  arbitrary  periodic  sequence  of  numbers  of  period  k 


(3.  4) 


It  is  now  easy  to  determine  the  monospline 

1  S(  z)  *  8 


K(z)  =  -±-  z  -  S(z)  , 


mi 


m,  k 


having  least  L^-norm,  or 


2  V 


ie. 


(3.5)  IIkII  =  (~  f  |K(  z)  J2  de)  2  =  minimum  ,  (z  =  e*“) 

2  \2tt  q  / 

Indeed,  by  Lemma  4  we  may  write 

M  z)  -  ~  zm  -  Y  b  .--b  .  zVt1 

m!  2m  y  v  v-m+1  v 

/ 1  1  ,  .  \  m  I  V  w  u  v 

Uil  2^i  m  1  my  2*1  l-1  v  v-m+1  v 
v  '  v#m 

and  using  Parseval’s  relation  we  obtain 


-h- 
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~  /  |k(  Z)  I2  d0  =  |-*T  -  ~~b  ...bn  I2  +  ^4-  E 

2"  ■>  1  I  |m!  2ni  m  1  'ml  ,2  u 


(3.6) 


4"  v=m(k) 
v#m 


v-m+1 


+  E  I’IjI2  E  |b„*  •  *  bv-m+l|2  • 

4tt  i=0  1  1  vs|(k)  1  v  m+i| 


Observe  that  by  (2.6) 

E  Iv  •  •  bv-m+l|2  >  0  for  all  integer  i  . 
vsi(k)  1  m+il 

It  follows  that  in  minimizing  the  right  side  of  (3.6)  we  must  have 


(3.7) 


=  0  if  I  =  0, 1, . . .  f  m-1,  m+1, . .. ,  k-1 


Therefore  the  search  reduces  to  the  much  narrower  class  of  monosplines  of  the 


(3.8) 


vsm(k) 


In  fact,  by  Lemma  3,  we  see  that  (3.8)  describes  the  most  general  mono¬ 
spline  K(  z)  satisfying  the  relation 


(3.  9) 
Setting 


K(  z  u)  =  w  K(z) 


we  may  write  (3.  8)  as 


(3.  10) 


ks+m 


(ks+l)(ks+z). .  .  (ks+m) 


and  (3.6)  becomes 


(3.11)  (m!  >2  (HkH  )2  =  ll-X  |  2  +  |\|2  E  - , - 7  , 

2  s*0  (ks+l)2(ks+2)^...(ks+m)2 

and  we  are  to  determine  the  value  of  X  that  minimizes  the  right  side.  The 
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inequalities  1 1-X  |  >  |l-Re  X  |  ,  |x|  >  | Re  V  |  show  that  this  value  of  X 


must  be  real.  Writing 

(3.12)  c=£  - - 

s*0  (ks+1 )  .. .  (ks+m) 

we  see  from  the  identity 

(1-X)2  +  cX2  =  (lt-c)(X-j^)2  +  >  (c>°>  » 

that  the  least  value  of  the  left-hand  side  is  =  c/(l+c)  and  that  it  is  assumed 
if  X  =  X+  =  i/(i +C)  .  By  (3.12)  we  see  that  the  right  side  of  (3.11)  is  least 
if  X  equals 

_ (m!)2 _ 

(ks+1)2  (ks+2)2  .  . .  (ks+m)2 

The  least  value  of  the  right  side  of  (3. 11)  we  now  find  to  be  equal  c/(l  +  c)  = 

(C-i'1.'1-1,  • 

We  have  therefore  established 

Theorem^.  The  unique  monospline  (3.  4)  of  least  L.-norm  is 

(3.14)  K*(  z)  =  *“T  zm  -  X+  S  (  z) 

*  m!  *  m 

where 

„  v  1  _ks+m 

(3-15)  Sm(Z)  "  It  (ks  +  l)(ks+2)  .  .  .  (ks  +  m)  2  ’ 

S  =  -oo 

while  X#  is  given  by  (3.13).  The  value  of  the  least  L.-norm  is 

(3.i6)  IM2  =  ■ 

4.  The_bes^2ii|^£Siiii™l2m«i«  (1  •  9)  .  We  obtain  the  best  quadrature 
formula  (1.9),  by  the  definition  of  the  term  "best"  introduced  by  A.  Sard  (  2], 
if  we  choose  in  (1 . 9)  K(  z)  =  K*  ( z)  ,  and  there  remains  to  determine  the 


(3.13) 


00 

£ 

is  =-00 
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coefficients  according  to  U.10).  However,  we  have  already  noticed 


that  K+(z)  satisfies  the  relation  (3.9)  which  by  differentiations  yields 


K*v\z  oj)  =  wm  VK+(z)  and  in  particular 


KT  "(zu)  =  «K'"*  "(z) 


rW-Ui 


By  (1.10)  we  see  that 

(4.1)  C  *Jco 

and  it  therefore  suffices  to  determine  the  jump  AQ  =  -  CQ  of  K^m  ^(z)  at 
z  =  1  .  It  is  easily  obtained  as  follows. 

Differentiating  the  Fourier  series  (3.  14)  m-1  times  we  get 

(4.2)  Kim-!,U>  =  z-X,£j-i-zkStl  • 

^  *  g  kS+l 

Its  sum  is,  as  we  know,  a  step  function  with  discontinuities  at  the  <J  .  To 
determine  its  Jump  AQ  at  z  =  1  ,  we  consider  the  function  <p(  0)  of  period 
2"  such  that 

(4.3)  ?(0)  =  i(ir-G)  if  0  <  0  <  2tt  . 


Defining  also  ?(0)  =  0  ,  we  have  the  familiar  expansion 

(4.4)  „(0)  =  f  sinw0  =  X  £  ieis0  for  all  real  0, 

1  v  s#0  s 

from  which  we  derive  that 


1  ....  10  1  v  1 

k  ^(k0)e  =  IT  2.  k7z 

s*0 


1  ks+1 
ks  Z 


( z  =  e 


i0 


Using  (4.  2)  we  find  the  expansion 

(4.  5)  K^m_1)(eie)  +  21  V*  I^(k0)ei0  =  (l-\„)  z  +  \  £ 

s*0 


which  is  seen  to  converge  absolutely  and  uniformly  on  U  . 


) 


I  ks+1 
ks(ksfl)  2 

We  conclude  that 
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the  left  side  belongs  to  C(  U)  .  The  jump  of  the  function  9)  being  =  w  , 
at  0  =  0,  by  (4.  3),  we  conclude  that  the  jump  2iX.%k  *tt  of  the  second 
term  of  the  left  side  of  (4.  5)  cancels  the  jump  AQ  of  the  first  term.  Therefore 

(4.6)  =  • 

Using  the  relations  (4.1)  we  obtain 

Theorem  2.  Let  1  s  m  <  k  .  Among  all  quadrature  formulae  of  the  form 
(1)  and  having  the  property  (1.  1),  the  best  is  the  Quadrature  formula 

(4.7)  /f(z)uz  =  X+2ffik_1  £  Jf(J)  +  (-l)"1  fKJz)fim)iz)dz 

z  J=0  U 

where  K+(z)  is  the  monospline  of  Theorem  1 .  while  the  positive  fraction  X. 
is  defined  by  (3.13)  . 

Theorem  2  does  not  state  or  imply  that  the  familiar  average  (or  Riemann 

sum) 

,k-i 

(4.8)  A=2Trik  £Jf{wJ ) 

0 

is  not  as  good  an  approximation  of  the  integral  as  the  modified  approximation 
\+A  appearing  in  (4.  7)  .  However,  Theorem  2  does  state  that  \+A  is  the 
best  approximation  in  the  following  sense:  For  all  functions  f(z)  such  that 

(4.9)  /  |f(m)U)I2de  <  H2  , 

0 

where  H  is  a  fixed  constant,  we  have  the  inequality 

(4.10)  |/f(z)dz  -  \%A|  S  J^-H  ,  J*  =  -VUMl-**))'  • 

U  m* 

For  any  other  inequality 
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k-1 

|  /  f (  z)dz  -  V  C  f(uJ)|  S  J  -  H 
U  0  3 
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valid  for  the  class  (4.  9) ,  we  must  have  J  >  J*  . 


®  i 

m,  k 


We  wish  to 


determine  the  quadrature  formula 


k-1 


(5.1) 

such  that 


/ f(  z)dz  =  Yj  C  f(  J)  +  ?f 


U 


(5. 2)  Rf  =  0  if  f  *  S  . 

m  y  K 

One  way  to  construct  such  a  formula  is  as  follows. 

S  ,  has  the  interpolatory  property  at  the  knots  J 
m,  < 

we  must 


We  assume  that  the  class 
.  By  [  4,  Theorem  1  ] 


(5.  3)  exclude  the  case  that  k  is  even  and  m  is  odd. 

We  are  then  assured  [4]  of  the  existence  of  a  fundamental  function 
Lm<z>'gm,*  such  that 

^  1  .  . 

(5.4)  f(z)  =  Y  f(«)L  »  v  • 

Integrating  both  sides  of  this  identity  along  U  we  see  that  (5.1)  has  the 
property  (5.  2)  if  we  set 

(5.  5)  C  =  /  L  (zu'1  )dz  . 

i  ^  m 

Again  we  assume  (5.  3)  .  A  second  way  of  constructing  this  quadrature 
formula  is  as  follows.  Let  SQ(t)  be  the  unique  spline  function  such  that 
the  monospline 

(5.6)  K0(z)  =  Jy  -  S0(z) 
has  the  property 

(5.7)  Kq(J)=0  (J  =  0,  ...,k-l)  . 
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The  existence  and  unicity  of  Sq(z)  follow  from  [4,  Theorem  1].  We  shall 
now  establish  several  properties  of  KQ(z)  which  will  show  that  this  mono¬ 
spline  will  produce  the  quadrature  formula  we  are  looking  for. 

I.  The  monos  pi  ine  (5.6)  satisfies 

(5.8)  Kq(zw)  =  J"K0(z)  . 

Proof:  Clearly 

-m„  .  ,  1  m  -m  _  .  . 

w  KQ  (  z  oj)  =  —  z  -  <*>  bQ(Zw) 

is  a  monospline  satisfying  the  relations  (5.7)  .  The  unicity  of  KQ(z)  impli  >" 
the  relation  (5.8)  . 

II.  The  quadrature  formula 

k-i 

(5.9)  / f(z)dz  =  £  C  nJ)  +  RQf  , 

U  0  3 

where  -  C°  =  the  lump  of  id™  ^(z)  at  z  =  J  .  is  exact  If  ft  8  ,  . 

-  j  - 1 - K -  O'  —  ’  -  m,  k 

Proof:  By  Lemma  1  we  may  write 

*  R0f  =  /  K  (z)f(m)(z)dz  =  f  K  (z)df(m_1)(z)  If  f€Cm(U>  . 

U  U 

However,  the  last  form 

(5.10)  *R  f  =  f  K  (z)df(rn_1)(z) 

U 

is  also  applicable  (see  proof  of  Lemma  1)  if  we  only  assume  f^m  ^(z)  to 

be  absolutely  continuous  and  f^m  ^(z)  of  bounded  variation  on  U  .  We 

may  therefore  assume  in  (5.9),  (5.10)  that  f(z)*S  ,  and,  therefore,  that 

m,  k 

f^m  ^(z)  is  a  step  function.  If  we  denote  by  6^  its  Jump  at  z  =  ,  then 

(5.7)  and  (5.10)  show  that 
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III.  C  =  cj  (J  =0,...,k-l)  . 

Proof:  Letting  f(z)  =  L  (zu  in  (5.9),  then  (5.5)  shows  that 


0  f  —  j  rv 

c  =  /  L  ( zco  J)d2  =  C  . 
j  m  j 


The  explicit  construction  of  the  quadrature  formula  (5.  9)  now  presents 
no  difficulties.  We  conclude  from  the  relation  (5.8)  that  SQ(z)  in  an 


m-flower  and  Lemma  3  implies  that  KQ(z)  may  be  written  as 


(5.11) 


v*>  *  ir zm  -  t 


S--00 


(ks+-l)(ks+2).  .  .  (ks+m) 


ks+m 


for  an  appropriate  value  of  the  constant  \Q  .  The  defining  properties  (5.7), 


in  particular  KQ ( 1 )  =  0  show  that 


(5.  12) 


s=-00 


(ks4-l)(ks+-2).  .  .  (ks+m) 


This  requires  that  the  sum  of  the  series  in  (5.12)  should  not  vanish.  By 
[4,  Lemma  4]  this  series  can  vanish  only  if  m  =  k-1  and  m  is  odd.  But 
this  would  imply  that  k  is  even  and  m  is  odd,  a  situation  which  can  not  occur 
because  of  our  assumption  (5.  3)  . 

We  summarize  our  results  in 
Theorem  3.  If  we  exclude  the  case  that 


(5.13) 


m  =  k-1  and  m  is  odd  , 


-i  v1..),,  i 


(5.14)  f  f(  z)dz  =  X.  2irik_l  £  Jf(J)  +  (-l)m  f  K.(  z)  f(m)(  z)dz 
U  0  U 

is  a  quadrature  formula  that  is  exact  if  f(z)«8  where  K  (z)  and  \n 

rrij  )c  0  0 

are  described  by  (5.11)  and  (5.12)  .  If  we  replace  (5.13)  by  the  more  stringent 


assumption  (5.3)  then  the  quadrature  formula  (5.14)  may  also  be  obtained  bv 


#1002 


-13- 


integration  from  the  spline  interpolation  formula  (5.4) 


The  value  CQ  =  2iriXQ/k  used  in  writing  the  formula  (5.14)  follow 
from  (4. 6)  on  replacing  X#  by  XQ  . 

6.  The  polynomial  components  of  a  certain  spline  function.  The  spline 


function 

(6.1) 


oo 


SJz)  =  2 


1 


ks+m 


(ks+l)(ks+2). . .  (ks+m) 


S=-oo 

appears  in  the  representation  (3.14),  (3.15)  of  K+(z)  as  well  as  in  (5.1)) 
which  describes  KQ(z)  .  By  Lemma  3  it  may  also  be  characterized  as  an 
m-flower  of  degree  m-1  .  Here  we  propose  to  determine  explicitly  its  in¬ 
dividual  polynomial  components.  Without  loss  of  generality  we  may  assume 
that  m  =  k-1  and  determine  the  components  of 


(6.2) 


<W>  *  l 


1 


ks+k-1 


-00 


(ks+1). .  .(ks+k-1) 


S  (z)  =S<k:m-1,(z)  ,  (lSmSk-1)  . 
m  k-l 


it  being  clear  by  differentiations  that 

(6.3) 

The  functional  equation 

(6.4) 

shows  that  it  suffices  to  determine  the  polynomial  component  in  the  arc  (1,  to) 
Setting 


S^_j(z  w)  =  u  S^_j(  z) 


(6.5) 


^  =  exp(nl/k) 


**1  Jc 

we  turn  Sk  ^(z)  by  ^  to  obtain,  since  4»  =  -1  ,  the  new  spline  function 


(6.6) 


00 

s(z)  =+S<*+)  =  l  .f-U 


s+1 


ks+k-1 


-to  (ks+1). ..  (ks+k-1) 


having  its  knots  in  the  new  locations  .  Let 
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k-2 


k-3 


(6.7) 


P0(z) 


- -  to.  — - -  +  +  a,  _ 

0  (k-2) 3  1  (k-3) I  k-2 


be  the  component  of  s(z)  in  the  arc  (ijT1,  iJj)  .  The  Laurent  series  (6.6) 
having  real  coefficients  it  is  clear  that  PU)  =  P(z)  if  z  is  on  the  arc 


( »|/)  and  therefore  also  for  all  z  .  This  shows  that  PQ(z)  is  a  real 
polynomial.  Also  s(z)  satisfies  the  relation  s(zo)  =w  *s(z)  whence 


(6.8) 


s(z)  =  CJ  ^s(zu)  ^  ) 


If  z  is  ir.  the  arc  (^,  i{ju>)  ,  then  (6.8)  and  (6.7)  show  that  the  com¬ 


ponent  of  s(z)  on  this  arc  is 

k-2 


k-3 


(6.9)  P,U>  =*„<■>  (tijj  to,-2  , • 

k-3 

At  this  point  we  use  that  s(z)€  C  (U)  and  obtain  at  z  =  ^  the  k-2 
equations 

(6.10)  PqV,)<4‘)  =  Pj  V')(4/>  (v  =  0,  ...,k-3)  . 

Using  (6.7)  and  (6.9)  we  obtain  the  system 

o0(u-l)  +  ^(w2-!)  7  +  •  •  •  +  av(uVH-l)  =0  (V  =1,  .  .  • ,  k-2) 

v  -4*1 

Dividing  by  w  -  1  and  using  the  relations 


a  , 
w  -i 


eg  —ct 

J _ 


_1 _  sin( ira/ k) 


p  .  ’  (p-cr)!  '  (p-a)l  ,P  -(3  ~  (p-ff)J  sin(TTp/k) 
w  -i  ~y 

we  obtain  the  system 

0  =  orQ  sin(ir/k)  f-  aj  sin(2"/k) 

0  =  ~  aQ  sin(rr/k)  +  sin(2n/k)  +  sin(3n/k) 


0  *  (k^TT  “o  sln(’/k)  *  (k^5)T  "i  * 


+  a-k2sin((k-l)ii/k) 
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which  in  terms  of  the  new  unknowns 


(6.11) 

becomes 


v  v  k 


0  =  IT  P0  +  P1 
0  =  IT  +  IT  pi  +  p2 


(v  =  0,  ..  .,k-2) 


0  = 


(k-2)l  P0  +  {k— 3)  1  P1  +  '  "  +  1)  Pk-3  +  Pk-2  ' 


This  being  the  system  giving  the  coefficients  of  the  rower  series  reciprocal 

x  v  V  v 

to  e  =  IjX  /v\  we  obtain  a  solution  =  (-1)  /v\  and  by  (6.11) 

(6.12)  %  =  (-i)/G!sin(i~“)  ,  (v  =  0, ..  .,k-2)  . 


Ftom  (6.6)  and  (6.7)  we  therefore  obtain  the  identity 


(6.13)  £  - 

'  i-i  /U  „  1  1  \ 


ks+k-1 


(ks+1). . .  (ks+k-1) 


1  =  c  .k£Vi>f;2) 

v=0 


k-v-2/ .  (v+l)n 

z  /sin  *-  -y1- 


-l 


valid  in  the  arc  (<\i  ,  40  ,  while  c  is  a  constant  yet  to  be  determined. 

To  obtain  c  we  set  z  =  1  to  obtain 

f.  O^^r,  ■ • -T< 

However,  the  left  side  is  easily  evaluated  directly  from  the  partial  fraction 
decomposition 
_ 1 


(6.15) 


_ _  1  f—k _ _1_+  f  (-i)k"2__l_\ 

:-l)  "  (k-2)  1  Hks+1)  k  1  J  ks+2  ‘  ks+k-lj 


(ks+1). . .  (ks+k- 

From  the  classical  expansion  £  (-l)s/(ii+x)  -  w/sinwx  we  obtain 


y  irif 

u  ks+v+1 


k  sin 


(vfl)ff 
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If  we  multiply  (6.15)  by  (-1)  and  sum  overall  s  we  have 


y _ Lzu! _ 

(ks+1). . .  (ks+k 


__  =  _jl _ kfZ(-nv(k~z) 

-1)  (k-2)  1  k  £  1  L}  {  v  J/Sln  k 


Comparing  with  (6.14)  we  find  that 


(6.16) 


(k-2)  lk  * 


Using  the  relation  (6.  3)  we  can  also  determine  the  polynomial  com¬ 


ponent  of  (6.1)  within  the  arc  (1,  to)  .  Indeed 


ks+m 


Sm( z^)  -  S  (ks+1).".1!  (ks+m)  2 
s 


and  differentiating  (6.13)  k-m-1  times  we  obtain  that 


*1  .. 


if  z  is  confined  to  the  arc  (4»  ,  «{*)  . 

As  an  example  let  us  determine  K+(z)  using  the  formulae  (3.14), 
(3.15)  for  k  =  5  and  m  =  2  .  The  rapidly  convergent  series  (3.13)  shows 
that  =  .969690  ,  and  (6.17)  furnishes  the  expression 


i(f  ZK.(zi|»)  =  \zZ  -  X 


*  5  sin( 


w/ 5)  (j 


-  * 
2  cos  — 


in  the  arc  (4»  40  •  Graphing  the  image  of  the  arc  (4>  *,40  by  the  quadratic 

2 

polynomial  on  the  right  side  and  turning  it  by  4>  =  w  (or  +72°)  we  obtain 
the  image  of  the  arc  (1,  to)  by  w  =  K+(z)  .  The  functional  equation 

z  w)  =  «ZK^(z) 

allows  to  complete  the  image  of  the  entire  circle  U  which  is  shown  in  Fig.  1 . 
The  five  corners  of  the  curve  are  the  images  of  the  a/  ,  Fig.  1  showing  in 


parentheses  the  corresponding  .  Moreover 
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